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Abstract 

We study the quantum evolution of black holes immersed in a 
de Sitter background space. For black holes whose size is comparable 
to that of the cosmological horizon, this process differs significantly 
from the evaporation of asymptotically flat black holes. Our model 
includes the one-loop effective action in the s-wave and large N ap- 
proximation. Black holes of the maximal mass are in equilibrium. 
Unexpectedly, we find that nearly maximal quantum Schwarzschild- 
de Sitter black holes anti-evaporate. However, there is a different per- 
turbative mode that leads to evaporation. We show that this mode 
will always be excited when a pair of cosmological holes nucleates. 
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1 Introduction 



Of the effects expected of a quantum theory of gravity, black hole radiance |l| 
plays a particularly significant role. So far, however, mostly asymptotically 
flat black holes have been considered. In this work, we investigate a qualita- 
tively different black hole spacetime, in which the black hole is in a radiative 
equilibrium with a cosmological horizon. 

The evaporation of black holes has been studied using two-dimensional 
toy models, in which one-loop quantum effects were included 0, ^, |]. We 
have recently shown how to implement quantum effects in a more realistic 
class of two-dimensional models, which includes the important case of di- 
mensionally reduced general relativity 0. The result we obtained for the 
trace anomaly of a dilaton-coupled scalar field will be used here to study the 
evaporation of cosmological black holes. 

We shall consider the Schwarzschild-de Sitter family of black holes. The 
size of these black holes varies between zero and the size of the cosmological 
horizon. If the black hole is much smaller than the cosmological horizon, the 
effect of the radiation coming from the cosmological horizon is negligible, and 
one would expect the evaporation to be similar to that of Schwarzschild black 
holes. Therefore we shall not be interested in this case. Instead, we wish to in- 
vestigate the quantum evolution of nearly degenerate Schwarzschild-de Sitter 
black holes. The degenerate solution, in which the black hole has the maxi- 
mum size, is called the Nariai solution 0. In this solution the two horizons 
have the same size, and the same temperature. Therefore they will be in 
thermal equilibrium. Intuitively, one would expect any slight perturbation of 
the geometry to cause the black hole to become hotter than the background. 
Thus, one may suspect the thermal equilibrium of the Nariai solution to be 
unstable. The initial stages of such a run-away would be an interesting and 
novel quantum gravitational effect quite different from the evaporation of an 
asymptotically fiat black hole. In this paper we will investigate whether, and 
how, an instability develops in a two-dimensional model derived from four- 
dimensional general relativity. We include quantum effects at the one-loop 
level. 

The paper is structured as follows: In Sec. |^ we review the Schwarzschild- 
de Sitter solutions and the Nariai limit. We discuss the qualitative expec- 
tations for the evaporation of degenerate black holes, which motivate our 
one-loop study. The two-dimensional model corresponding to this physi- 
cal situation is presented in Sec. |^, and the equations of motion are de- 
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rived. In Sec. ^ the stability of the quantum Nariai solution under different 
types of perturbations is investigated. We find, quite unexpectedly, that the 
Schwarzschild-de Sitter solution is stable, but we also identify an unstable 
mode. Finally, the no-boundary condition is applied in Sec. |^ to study the 
stability of spontaneously nucleated cosmological black holes. 



2 Cosmological Black Holes 
2.1 Metric 

The neutral, static, spherically symmetric solutions of the Einstein equa- 
tion with a cosmological constant A are given by the Schwarzschild-de Sitter 
metric 

ds'^ = -V{r)dt^ + V{r)-^dr^ + r^dn\ (2.1) 

where 

V{r) = 1 - ^ - ^,2. (2.2) 

dVt^ is the metric on a unit two-sphere and is a mass parameter. For 
< /i < |A~^/^, V has two positive roots and rb, corresponding to the 
cosmological and the black hole horizons, respectively. The limit where /i — >^ 
corresponds to the de Sitter solution. In the limit /i — >■ |A~^/^ the size of the 
black hole horizon approaches the size of the cosmological horizon, and the 
above coordinates become inappropriate, since V{r) — > between the two 
horizons. Following Ginsparg and Perry |^, we write 

g^u^A = 1 - 3e^ < e < 1. (2.3) 

Then the degenerate case corresponds to e ^ 0. We define new time and 
radial coordinates ip and x by 



1 , 1 
r = — i=ip 



A VA 



1 2 

1 - ecosx - -e 
o 



(2.4) 



In these coordinates the black hole horizon corresponds to x = and the 
cosmological horizon to x = vr. The new metric obtained from the transfor- 
mations is, to first order in e, 

ds^ = (^1 + ^ecosx) sin\ (i^/-^ + ^ (^1 - ^ecosx) (ix^ (2.5) 
-t- ^ (1 - 2ecosx)rffi2- 
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This metric describes Schwarzschild-de Sitter solutions of nearly maximal 
black hole size. 

In these coordinates the topology of the spacelike sections of Schwarzschild- 
de Sitter becomes manifest: x S^. In general, the radius, r, of the 
two-spheres varies along the 5*^ coordinate, with the minimal (maximal) 
two-sphere corresponding to the black hole (cosmological) horizon. In the 
degenerate case, the two-spheres all have the same radius. 

2.2 Thermodynamics 

The surface gravities of the two horizons are given by 




where the upper (lower) sign is for the cosmological (black hole) horizon. In 
the degenerate case, the two horizons have the same surface gravity, and, 
since T = h/2ti, the same temperature. They are in thermal equilibrium; 
one could say that the black hole loses as much energy due to evaporation as 
it gains due to the incoming radiation from the cosmological horizon. Away 
from the thermal equilibrium, for nearly degenerate Schwarzschild-de Sitter 
black holes, one could make the simplifying assumption that the horizons still 
radiate thermally, with temperatures proportional to their surface gravities. 
This would lead one to expect an instability: By Eq. (|2.6|) , the black hole will 
be hotter than the cosmological horizon, and will therefore suffer a net loss of 
radiation energy. To investigate this suspected instability, a two-dimensional 
model is constructed below, in which one-loop terms are included. 

3 Two-dimensional Model 

The four-dimensional Lorentzian Einstein-Hilbert action with a cosmological 
constant is 

lovr J 

where R}^ and g^^ are the four- dimensional Ricci scalar and metric determi- 
nant, and the fi are scalar fields which will carry the quantum radiation. 



N 



/?^v_2A-^^(Viv^, 

^ i=i 



(3.1) 
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We shall consider only spherically symmetric fields and quantum fiuctu- 
ations. Thus, we make a spherically symmetric metric ansatz, 



ds^ = e'p (-de + dx") + e-^^dn^ 



(3.2) 



where the remaining two-dimensional metric has been written in conformal 
gauge; x is the coordinate on the one-sphere and has the period 27r. Now the 
spherical coordinates can be integrated out, and the action is reduced to 



S=— (d^xi-gf'h 

1671 J ^ 



-20 



N 



R + 2(V0)2 + 2e^^ -2A- ^(V/^; 



i=l 



(3.3) 



where the gravitational coupling has been rescaled into the standard form. 
Note that the scalar fields have acquired an exponential coupling to the 
dilaton in the dimensional reduction. In order to take quantum effects into 
account, we will find the classical solutions to the action 5* + W*. W* is 
the scale-dependent part of the one-loop effective action for dilaton coupled 
scalars, which we derived in a recent paper p: 



W* 



A8n 



dM 



a/2 



-R-R - 6(V0)2-i? - 20i? 



(3.4) 



The (V0)^ term will be neglected; we justify this neglect at an appropriate 
place below. 

Following Hayward |P, we render this action local by introducing an 
independent scalar field Z which mimics the trace anomaly. The / fields 
have the classical solution fi = and can be integrated out. Thus we obtain 
the action 



S 



IGtt 



jd\{-gY'' 



R 



K 



-{VZY + 2 



2e-^^ (V0)^ - 2e-^^A 



where 



2N 



K = 



(3.5) 



(3.6) 



There is some debate about the coefficient of the (pR term in the effective 
action. Our result corresponds to the choice w = 2; the RST coefficient |§] 
corresponds to w = 1, and the result of Nojiri and Odintsov |10| can be 
represented by choosing w = —6. In Ref. [0, probably erroneously, w = —1 
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was chosen. We take the large N hmit, in which the quantum fluctuations 
of the metric are dominated by the quantum fluctuations of the scalars; 
thus, K ^ 1. In addition, for quantum corrections to be small we assume 
that b = kA <^ 1. To first order in b, we shall find that the behaviour of the 
system is independent of w. 

For compactness of notation, we denote differentiation with respect to t 
(x) by an overdot (a prime). Further, we define for any functions / and g: 

dfdg^-fg + f'g\ d'g^-g + g", (3.7) 

and 

5f5g^fg + fg\ 6'g^g + g". (3.8) 

Variation with respect to p, cf) and Z leads to the following equations of 
motion: 

1 _ !^e2'^) 9^0 + 2{d<pf + '^e^'t'd'^Z + e2^+20 ^^^~24> _ i) = Q; (3.9) 

d^Z-2d^p = {). (3.11) 
There are two equations of constraint: 

1 _ (^2^ _ _ (50)2 ^ ^^24> U^zf + 25^ Z - ASZ5p 

4 / V / g L 

(3.12) 

1 _ !^e2*) (0' - p0' - p'0) - 00' = ^e^^ [ZZ' + 2Z'-2 [pZ' + p' z) . 

(3.13) 

From Eq. ( |3.11| ), it follows that 

Z = 2p + i], (3.14) 

where i] satisfies 

d'^rj = 0. (3.15) 

The remaining freedom in rj can be used to satisfy the constraint equations 
for any choice of p, p, and on an initial spacelike section. This can 
be seen most easily by decomposing the fields and the constraint equations 
into Fourier modes on the initial S^. By solving for the second term on 
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the right hand side of Eq. ( |3.12| ), and by using Eqs. ( ^.14| ) and ( ^.15| ), the 
first constraint yields one algebraic equation for each Fourier coefficient of 
rj. Similarly, the second constraint yields one algebraic equation for the 
time derivative of each Fourier coefficient of 1]. If the initial slice was non- 
compact, this argument would suffice. Here it must be verified, however, 
that rj and fj will have a period of 2tt. The problem reduces to the question 
whether the two constant mode constraint equations can be satisfied. Indeed, 
while for each oscillatory mode of rj, there are two degrees of freedom (the 
Fourier coefficient and its time derivative), the second time derivative of the 
constant mode coefficient, f)o, must vanish by Eq. ( |3.15| ). Thus there is only 
one degree of freedom, -fiQ, for the two constant mode equations. However, 
since we have introduced no odd modes (i.e., modes of the form sinkx) in 
the perturbation of (j), none of the fields will contain any odd modes. Since 
each term in Eq. ( |3.13| ) contains exactly one spatial derivative, each term 
will be odd. Therefore all even mode components of the second constraint 
vanish identically. In particular the constant mode component will thus be 
automatically satisfied. Then the freedom in 7)0 can be used to satisfy the 
constant mode component of the remaining constraint, Eq. ( |3.12|) , through 
the first ^ term on the right hand side. 



4 Perturbative Stability 
4.1 Perturbation Ansatz 

With the model developed above we can describe the quantum behaviour of a 
cosmological black hole of the maximal mass under perturbations. The Nariai 
solution is still characterised by the constancy of the two-sphere radius, e~'^. 
Because of quantum corrections, this radius will no longer be exactly A^^/^. 
Instead, the solution is given by 

^"' = ^^,^ ^'' = ^2, (4.1) 
Ai cos^t 

where 

1 1 



Ai 8A 



4 - {w + 2)b + a/16 - 8(w-2)6 + (^ + 2)252 



(4.2) 



"'^Note that rjo can thus be purely imaginary, as indeed it will be for the Nariai solution, 
signaling negative energy density of the quantum field. 
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A, ' 



2WK 

Expanding to first order in 6, one obtains: 



4 + {w + 2)h - ^lQ-S{w-2)h+{w + 2)W 



(4.3) 




(4.4) 

Let us now perturb this solution so that the two-sphere radius, e~'^, varies 
shghtly along the one-sphere coordinate, x: 

e^^ = A2 [1 + 2ta{t) cos x] , (4.6) 

where we take e -C 1. We will call a the metric perturbation. A similar 
perturbation could be introduced for e^'', but it does not enter the equation 
of motion for a at first order in e. This equation is obtained by eliminating 
d^Z and d"^ p from Eq. (|3.9|) using Eqs. (p.ll|) and (|3.1CI|) , and inserting the 
above perturbation ansatz. This yields 

^ = ^ - 1. (4.7) 

where 



2^16 - 8{w-2)b + (w + 2)262 



A — wb 



(4.8) 



To first order in 6, one finds that 



2 + 6, (4.9) 



which means that w, and therefore the (pR term in the effective action, play 
no role in the horizon dynamics at this level of approximation. This is also the 
right place to discuss why the term ^/—g (VcpY^R in the effective action can 
be neglected. In conformal coordinates this term is proportional to {dipYp. 
Thus, in the p-equation of motion, Eq. ( |3.9| ), it will lead to a (90)^ term, 
which is of second order in e and can be neglected. In the 0-equation of 
motion, Eq. (|3.10| ), it yields terms proportional to k that are of first order 



in e. They will enter the equation of motion for a via the ne^'^d'^Z term in 
Eq. ( p.lO| ). Thus they will be of second order in b and can be dropped. The 
neglect of the log /i^ term |0 can be justified in the same way. 
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4.2 Horizon Tracing 

In order to describe the evolution of the black hole, one must know where 
the horizon is located. The condition for a horizon is (V0)^ = 0. Eq. (4.6) 
yields 

d(f) d(f) 

^r- = ecrcosx, T^ = —easmx. (4.10) 
ot ox 

Therefore, the black hole and cosmological horizons are located at 



Xb(t) = arctan 



Xe(t) = 7r-Xb(t). (4.11) 

To first order in e, the size of the black hole horizon, rb, is given by 

rb(t)-2 = e^'^t*'^''^^] = As [1 + 2e6{t)] , (4.12) 
where we define the horizon perturbation 

(5 = cosxb = a ( 1 + — j . (4.13) 

We will focus on the early time evolution of the black hole horizon; the 
treatment of the cosmological horizon is completely equivalent. 

To obtain explicitly the evolution of the black hole horizon radius, rb(t), 
one must solve Eq. (p7| ) for cr(t), and use the result in Eq. ( [4. 13]) to evaluate 



Eq. ( 4.121) . If the horizon perturbation grows, the black hole is shrinking: 



this corresponds to evaporation. It will be shown below, however, that the 
behaviour of 5{t) depends on the initial conditions chosen for the metric 
perturbation, o"o and ao. 

4.3 Classical Evolution 

As a first check, one can examine the classical = 0. This has a = 2, 



and Eq. (4/7) can be solved exactly. From the constraint equations, Eq. ( 3.12 ) 
and ( p.l3| ), it follows that 

a = a taut. (4.14) 

Therefore the appropriate boundary condition at t = is (Tq = 0. The 
solution is 

a{t) = ^ (4.15) 
cost 
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Then Eq. yields 

6(t) = (Jq = const. 



(4.16) 



Since the quantum fields are switched off, no evaporation takes place; the 
horizon size remains that of the initial perturbation. This simply describes 
the case of a static Schwarzschild-de Sitter solution of nearly maximal mass, 
as given in Eq. (|2.6| ). 



4.4 Quantum Evolution 



When we turn on the quantum radiation (k > 0) the constraints no longer 
fix the initial conditions on the metric perturbation. There will thus be two 
linearly independent types of initial perturbation. The first is the one we 
were forced to choose in the classical case: cxo > 0, do = 0. It describes 
the spatial section of a quantum corrected Schwarzschild-de Sitter solution 
of nearly maximal mass. Thus one might expect the black hole to evaporate. 



For a > 2, Eq. (O) cannot be solved analytically. Since we are interested in 
the early stages of the evaporation process, however, it will suffice to solve 
for cr as a power series in t. Using Eq. ( |4.13| ) one finds that 



m 



(^0 



l--(a-l)(a-2)t2 + 0(t' 



1 



-bt^ 
2 



(4.17) 



The horizon perturbation shrinks from its initial value. Thus, the black hole 
size increases, and the black hole grows, at least initially, back towards the 
maximal radius. One could say that nearly maximal Schwarzschild-de Sitter 
black holes "anti-evaporate" . 

This is a surprising result, since intuitive thermodynamic arguments would 
have led to the opposite conclusion. The anti-evaporation can be understood 
in the following way. By specifying the metric perturbation, the radiation 
distribution of the Z field is implicitly fixed through the constraint equations, 
(|3.12| ) and ( p.l3|) . Our result shows that radiation is heading towards the 
black hole if the boundary condition ctq > 0, do = is chosen. 

Let us now turn to the second type of initial metric perturbation: ctq = 0, 
do > 0. Here the spatial geometry is unperturbed on the initial shce, but it 
is given a kind of "push" that corresponds to a perturbation in the radiation 
bath. Solving once again for a with these boundary conditions, and using 



10 



Eq. ( 4.13 ), one finds for small t: 

S{t) = cTot'^. (4.18) 

The horizon perturbation grows. This perturbation mode is unstable, and 
leads to evaporation. 

We have seen that the radiation equilibrium of a Nariai universe displays 
unusual and non-trivial stability properties. The evolution of the black hole 
horizon depends crucially on the tj^e of metric perturbation. Neverthe- 
less, one may ask the question whether a cosmological black hole will typi- 
cally evaporate or not. Cosmological black holes cannot come into existence 
through classical gravitational collapse, since they live in an exponentially 
expanding de Sitter background. The only natural way for them to appear is 
through the quantum process of pair creation 0. This pair creation process 
can also occur in an inflationary universe, because of its similarity to de Sit- 
ter space p, |lT], |12|. The nucleation of a Lorentzian black hole spacetime is 
described as the analytic continuation of an appropriate complex solution of 



the Einstein equations, which satisfies the no boundary condition |13|. We 
will show below that the no boundary condition selects a particular linear 
combination of the two types of initial metric perturbation, thus allowing us 
to determine the fate of the black hole. 



5 No Boundary Condition 

To obtain the unperturbed Euclidean Nariai solution in conformal gauge, 
one performs the analytic continuation t = ir in the Lorentzian solution, 



Eq. (|1]). This yields 

(ds^v^' = e^P (dr^ + dx") + e-^U^\ (5.1) 

and 

In four dimensions, this describes the product of two round two-spheres of 
slightly different radii, A^^ ^''^ and Ag^^^. The analytic continuation to a 
Lorentzian Nariai solution corresponds to a path in the r plane, first along 
the real r axis, from r = — ootor = 0, and then along the imaginary axis 
from t = Otot = tt/2. This can be visualised geometrically by cutting 
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the first two sphere in half, and joining to it a Lorentzian 1 + 1-dimensional 
de Sitter hyperboloid. Because the (r, x) sphere has its north (south) pole 
at r = oo (r = — oo), it is convenient to rescale time: 



sinw = — - — , (5.3) 
coshr 



or, equivalently, 



dr 

cosM = — tanhr, cot m = — sinh r, du = — . (5.4) 

coshr 

With the new time coordinate u, the solution takes the form 

(ds'^^y = 1- [du^ + siTi\dx^) + ^rffil (5.5) 

Now the south pole lies at u = 0, and the nucleation path runs io u = 7r/2, 
and then parallel to the imaginary axis (u = 7r/2 + iv) from f = to f = oo. 

The perturbation of e^*^, Eq. ( [4.6|) introduces the variable a, which satis- 
fies the Euclidean version of Eq. 



2 d da f 2 \ 

sin u -—r + sin u cos u— 1 — a sin m a = 0. (5.6) 

du^ du ^ ^ 

In addition, the nature of the Euclidean geometry enforces the boundary 
condition that the perturbation vanish at the south pole: 

a{u = 0) = 0. (5.7) 

Otherwise, e^"^ would not be single valued, because the coordinate x degen- 
erates at this point. This leaves a as the only degree of freedom in the 
boundary conditions at u = 0. 

It will be useful to define the parameter c by the relation c(c + 1) = a. 
The classical CclSG, Qj — 2, corresponds to c = 1; for small 6, they receive 
the quantum corrections a = 2 + 6 and c = 1 + 6/3. With the boundary 
condition, Eq. (|5.7| ), the equation of motion for a, Eq. (|5.6| ), can be solved 
exactly only for integer c (a = 2, 6, 12, 20,. . . ). The solution is of the form 

(^{u) = ^ Afc sin(c - 2fc)u, (5.8) 

0<fc<c/2 
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with constants Aj.. Even for non-integer c, however, this turns out to be a 
good approximation in the region < u < tt/2 of the {u,v) plane. Since we 
are interested in the case where 6 -C 1, the sum in Eq. ( |5.8| ) contains only 
one term, and we use the approximation^ 

a{u) ^ A sin cu. (5-9) 

It is instructive to consider the classical case first. (Physically, this is 
questionable, since the no boundary condition violates the constraints at 
second order in e.) For 6 = 0, the solution cr{u) = Asinu is exact. Along 
the Lorentzian line {u = tt/2 + iv), this solution becomes a{v) = A cosh v. 
By transforming back to the Lorentzian time variable t, one can check that 
this is the stable solution found in the previous section, with cxo = A, &o = 0. 
For real A, it is real everywhere along the nucleation path. Thus, when 
the quantum fields are turned off, the Euclidean formalism predicts that the 
unstable mode will not be excited. This is a welcome result, since there are 
no fields that could transport energy from one horizon to another. 

Once b is non-zero, however, it is easy to see that da/du no longer vanishes 
at the origin of Lorentzian time, u = tt/2. This indicates that the unstable 
mode, do 7^ 0, will be excited. Unfortunately, checking this is not entirely 
straightforward, because a is no longer real everywhere along the nucleation 
path. One must impose the condition that a and a be real at late Lorentzian 
times. We will first show that this can be achieved by a suitable complex 
choice of A. One can then calculate the horizon perturbation, 6, from the 
real late-time evolution of the metric perturbation, a, to demonstrate that 
evaporation takes place. 



^Treating Eq. ([5.6D perturbatively in b around a = 2 leads to untractable integrals. 



Fortunately the guessed approximation in Eq. (5.9) turns out to be rather accurate, espe- 
cially for late Lorentzian times v, which is the regime in which we claim our results to be 
valid. It is easy to check numerically that for sufficiently large v {v > 10), both the real 



and the imaginary part of Eq. (5.9) have a relative error 6/30 or less. The result for the 
phase of the prefactor, Eq. (5.13), has a relative error of less than 10~^, independently of 
b. Crucially, the exponential behaviour at late Lorentzian times is reproduced perfectly, 
as the ratio 

da/dv 
a 

using the approximation, agrees with the numerical result to machine accuracy. Therefore 



the relative error in Eq. (5.15) is the same as in Eq. (5.9); in both equations it is located 
practically entirely in the magnitude of the prefactor. — These statements hold for < 
b < I, which really is a wider interval than necessary. 
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From Eq. (|5.9|) one obtains the Lorentzian evolution of a, 



a{v) = Asmc(^^ + ivj (5.10) 

~ / CTT CTT \ 

= A [sin — cosh cv + i COS — siirn cv j . (5-11) 

For late Lorentzian times (i.e., large v), coshcf ~ sinhcf ~ e™/2, so the 
equation becomes 

(T(t;) ^ (ze-'^"/2^ e^^ (5.12) 
This can be rendered purely real by choosing the complex constant A to be 

i = A (-ie^""/2^ , (5.13) 

where A is real. 

Now we can return to the question whether the Euclidean boundary con- 
dition leads to evaporation. After transforming the time coordinate, the 
expression for the growth of the horizon perturbation, Eq. (|4.13| ), becomes 



6{v) = a 



1 + cosh V 



I da /di 



21 -1/2 



(5.14) 



The late time evolution is given by o{y) = 'fe'^". This yields, for large v, 

6{v) ^ -e™ (l + cV)"'^' ^ — exp(-t;). (5.15) 

This result confirms that pair created cosmological black holes will indeed 
evaporate. The magnitude of the horizon perturbation is proportional to the 
initial perturbation strength, A. The evaporation rate grows with kA. This 
agrees with intuitive expectations, since k measures the number of quantum 
fields that carry the radiation. 



6 Summary 

We have investigated the quantum stability of the Schwarzschild-de Sitter 
black holes of maximal mass, the Nariai solutions. From four-dimensional 
spherically symmetric general relativity with a cosmological constant and 
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minimally coupled scalar fields we obtained a two-dimensional model in 
which the scalars couple to the dilaton. The one-loop terms were included 
in the large limit, and the action was made local by introducing a field Z 
which mimics the trace anomaly. 

We found the quantum corrected Nariai solution and analysed its be- 
haviour under perturbations away from degeneracy. There are two possible 
ways of specifying the initial conditions for a perturbation on a Lorentzian 
spacelike section. The first possibility is that the displacement away from 
the Nariai solution is non-zero, but its time derivative vanishes. This per- 
turbation corresponds to nearly degenerate Schwarzschild-de Sitter space, 
and somewhat surprisingly, this perturbation is stable at least initially. The 
second possibility is a vanishing displacement and non-vanishing derivative. 
These initial conditions lead directly to evaporation. The different behaviour 
of these two types of perturbations can be explained by the fact that the 
initial radiation distribution is implicitly specified by the initial conditions, 
through the constraint equations. 

If neutral black holes nucleate spontaneously in pairs on a de Sitter back- 
ground, the initial data will be constrained by the no boundary condition: 
it selects a linear combination of the two types of perturbations. By finding 
appropriate complex compact instanton solutions we showed that this con- 
dition leads to black hole evaporation. Thus neutral primordial black holes 
are unstable. 
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